Possible high T c superconductivity mediated by antiferromagnetic spin fluctuations in 

systems with Fermi surface pockets 
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We propose that if there are two small pocket-like Fermi surfaces, and the spin susceptibility is 
pronounced around a wave vector Q that bridges the two pockets, the spin-singlet superconductivity 
mediated by spin fluctuations may have a high transition temperature. Using the fluctuation ex- 
change approximation, this idea is confirmed for the Hubbard on a lattice with alternating hopping 
integrals, for which T c is estimated to be almost an order of magnitude larger than those for systems 
with a large connected Fermi surface. 

PACS numbers: 74.20-z, 74.20Mn 
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Fascination for electronic mechanisms of superconduc- 
tivity has a long history. Given the large energy scale 
of the electrons itself, there has been an expectation 
from the early days that purely electronic mechanisms 
of superconductivity may lead to high transition temper- 
ature (T c ). The discovery of high T c superconductivity 
in the cupratescl has raised renewed interest in this pos- 
sibility. In fact, the repulsive Hubbard model, one of 
the simplest purely electronic models for the cuprates, 
has been studied extensively to investigate whether the 
model can account for T c ~ 100K. In particular, the fluc- 
tuation exchange (FLEX) studies have shown that the 
Hubbard model near half-filling on a two dimensional 
(2D) square lattice exhibits spin-singlet d-wave super- 
conductivity, mediated by antiferromagnetic spin fluctu- 
ations, having T c ~ O(0.01t)J3 Here t is the hopping 
integral in the Hubbard model, which is estimated to 
be ~ O(leV) for the cuprates, so that T c ~ 100 K, 
namely, a 'high T c '. Studies along this line have been 
performed on various lattices near half-filling, where T c 
of the (i-wave superconductivity has always turned out 
to be ~ O(O.Olt) on 2DQ"fLnd quasi ID latticesJHl and 
even lower on 3D lattices.113 

On the other hand, the present authors have recently 
studied the Hubbard model for various cases (i.e., com- 
binations of 2D or 3D lattice structures and band filling) 
having ferromagnetic spin fluctuations to find that spin- 
triplet superconductivity mediated by ferromagnetic spin 
fluctuations, if any, has very low T c in general.113 In fact, 
the possibility of finite T c has been suggested only in 
systems having two disconnected pieces of Fermi, surface 
lying point symmetrically about the T point .Ell There, 
superconductivity is enhanced because the nodal lines of 
the gap function do not intersect the Fermi surface. 

Conversely, we may say that a reason why the singlet 
d-wave superconductivity mediated by antiferromagnetic 
spin fluctuations has such 'low T c ' (~ O(0. 01t)) com- 
pared to the original electron energy scale t is because the 
gap function has nodes on the Fermi surface. Let us first 



see why this is the case. Superconductivity arises due 
to pair scattering processes near the Fermi surface me- 
diated by a pairing interaction. Within the BCS theory, 
contributions to superconductivity from pair scattering 
processes [k, -k] e FS -> [k', -k'] G FS (FS stands for 
Fermi surface) are summed up in the form 
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where V^(k. — k') is the pairing interaction and </>(k) is 
the superconducting gap function. In order to have su- 
perconductivity with an order parameter <f>, V c s has to 
be positive and large. When the pairing is mediated by 
spin fluctuations that are pronounced around a certain 
wave vector Q, the pairing interaction V^ 2 ' is positive and 
roughly proportional to the spin susceptibility, so only 
the pair scattering processes [k, — k] £ FS — > [k',— k'] £ 
FS that accompany a momentum transfer ~ Q and a 
sign change in <f> give significant positive contributions 
to Veff . This is the origin of the nodes in the d-wave gap 
function. The nodes intersect the Fermi surface for a con- 
nected Fermi surface, so that some pair scatterings give 
negative contributions to V c g as far as the spin suscepti- 
bility has a finite spread AQ (see Fig.|l](a)). In order to 
reduce negative contributions, the spin susceptibility has 
to have a sharp structure around Q, but in that case, pair 
scattering processes giving positive contributions will be 
strictly restricted to certain combinations of k and k', so 
in any case, V e f{ should be rather limited when the nodes 
of the gap function intersect the Fermi surface. 

Due to the above considerations, we propose that if (i) 
there are two small pocket-like Fermi surfaces A and B 
and (ii) the spin susceptibility is pronounced around a 
wave vector Q with a spread AQ, where Q bridges the 
Fermi surfaces A and B, and AQ is about the size of 
the pockets (see FigJ^(b)), the resulting spin-singlet su- 
perconductivity may have a T c much larger than 0.0 It. 
Namely, with different signs of the gap function assigned 
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FIG. 1. The d-wave gap function for a system with a 
large connected Fermi surface (a) or a gap function for sys- 
tems with two pocket-like Fermi surfaces A and B (b). The 
solid (dashed) curves represent positive (negative) values of 
the gap function on the Fermi surface. The solid (dashed) ar- 
rows represent spin-fluctuation-mediated pair scattering pro- 
cesses which give positive (negative) contribution to V c g. Q is 
the wave vector at which the spin susceptibility peaks, while 
AQ is the spread of the peak structure. 



on A and B (with no nodes on the Fermi surface), 
interpocket pair scattering processes V[k, — k] € A — > 
V[k',— k'] S B (and vice versa) all have significant pos- 
itive contribution to V e s, resulting in a high T c . This 
idea can be applied either to single band systems having 
disconnected pieces of Fermi surface or to systems with 
multiple Fermi surfaces corresponding to different bands. 

In order to provide an example in which such a sce- 
nario is realized, we start with a 2D rectangular lattice 
with different hopping integrals between the x and y di- 
rections (t X 2 — t x i = t x in Fig.@). For t x /t v > 1 the 
Fermi surface becomes open in the y direction (Fig.|3|,left 
panel). When the band filling is close to half- filling, i.e. 
n ~ 1, the Fermi surface intersects the lines k x = tv, 3tt/2 
for values of t x /t y not too largeJl3 The nesting vector Q 
in this case bridges the parts of the Fermi surface that 
lie inside and outside of the region tt/2 < k x < 3ir/2. If 
we now introduce an alternation in the hopping integrals 
in the x direction (t x \ ^ t X 2 in Figj^), the band will be 
folded at k x — n, 37r/2, resulting in two pocket-like Fermi 
surfaces, each in different bands (Fig.||,right panel). The 
important point here is that the nesting vector Q should 
not change, at least when t x \jt x % is not too far from 
unity, so that Q bridges the two Fermi surfaces. 
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FIG. 3. The Brillouin zone and the Fermi surface(s) for 
t x i = tx2 (right panel) and for t x \ ^ t x i (left panel). 



We have performed FLEX calculation for the Hub- 
bard model, H = -E<ij> =T ,lM c L c j<T 
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U^2n Uiiriii, op the above lattice. In the two band ver- 
sion of FLEXjj the Green's function G, the susceptibility 
X, the self-energy E, and the superconducting gap func- 
tion <{> all become 2x2 matrices, e.g., G/ m (k, is n ), where 
l,m specify the two sites in a unit cell. The orbital- 
indexed matrices for Green's function and the gap func- 
tions can be converted into band (denoted as a and 0)- 
indexed ones with a unitary transformation. As for the 
spin susceptibility, we diagonalize the 2x2 matrix Xsp 
and concentrate on the larger eigenvalue, denoted as x- 

The actual calculation proceeds as: (i) Dyson's equa- 
tion is solved to obtain the renormalized Green's function 
G(fc), where k is a shorthand for the wave vector k and 
the Matsubara frequency, ie n ; (ii) The effective interac- 
tion vM{q) is given as V^{q) = |y sp (g) + \V ch {q) - 
U 2 Xin(q)- Here, the effective interactions due to spin 
fluctuations (sp) and those due to charge fluctuations 
(ch) have the forms V sp = U 2 x sp and V ch = U 2 Xch, 
respectively, where the spin and the charge suscepti- 
bilities are x sp(ch )(g) = Xirr(g)[l - (+)L/xirr('7)] _1 in 
terms of the irreducible susceptibility matrix Xiniq) = 
— -^-^ fc G(fc + q)G(k) (A^:number of fc-point meshes); 
(iii) V"W then brings about the self-energy, S(/c) = 
jj J2 q G(k—q)V^ 1 ' (q), which is fed back to Dyson's equa- 
tion, and the self-consistent iterations are repeated until 
convergence is attained. We take 64 x 64 fc-point meshes 
and up to 4096 Matsubara frequencies. 

We determine T c as the temperature at which the 
eigenvalue A of the Eliashberg equation, 
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FIG. 2. The lattice considered in the present study. 



reaches unity. Here the pairing interaction for sin- 



glet pairing is given by V' 2 ' (q) 
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;V ch (q). 



We take t v as the unit of the energy (t y = 1). We 
present the results in the multiband scheme even for the 



single band case of t x i = t X 2- Throughout the study, the 
band filling (number of electrons / number of sites) and 
the on-site repulsion are fixed at n = 0.95 and U = 7, 
respectively. 

Let us now proceed to the calculation results. We 
first vary t x \ fixing t x % = 1.3 (solid arrow in the inset 
of Fig.|(a)). In Fig.|, we plot |G Q | 2 , \G \ 2 and X for 
the lowest Matsubara frequency for t x \ = 1.3 (a) and 
0.8 (b). The ridges of \G\ 2 represent the Fermi surfaces. 
For t x i = 1.3 (= t x2 ), this is a single band system, so 
the Fermi surfaces in the a and (3 bands are connected 
at around k = (-k,tt/2) and (tt, 3tt/2) to result in a sin- 
gle, large (open) Fermi surface. The spin susceptibility 
peaks at the nesting vector of the Fermi surface, namely 
at Q = (0, 7r) in the folded Brillouin zone (or Q — (tt, tt) 
in the unfolded Brillouin zone). The corresponding gap 
function for the lowest Matsubara frequency shown in 
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FIG. 4. Contour plots of |G a (k)r, |Ga(k)| 2 and X (q) for 
t xl =1.3(a) and 0.8(b) with t x2 = 1.3 and T = 0.14. The 
lighter the color, the larger the value. Note that \G a \ 2 in (b) 
is a contour plot of Fig. Me). 




FIG. 5. Plots of (p a (k) and 0/3 (k) for the same sets of pa- 
rameter values as in Figjl 
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FIG. 6. A(a) and T c (h) plotted as functions of t x i with 
T = 0.14 in (a) and t X 2 fixed at 1.3. T c is not obtained 
for t x i > 1.1, where the system comes very close to antiferro- 
magnetic ordering (U Xirr(Q) ~ 1) before A reaches unity upon 
lowering the temperature. The hatched region in (b) represent 
the range of T c for systems having a large connected Fermi 
surface. The solid(dashed) arrow in the inset of (a) shows the 
range of the parameter values taken in FigJ^(Fig.g|) . 



Fig.||(a) has a node along k x ~ n in order to satisfy the 
condition </>(k) = — </>(k + Q) on the Fermi surface. 

For t x \ = 0.8 by contrast, there is a small pocket-like 
Fermi surface in each band, namely, around k = (tt, tt) 
in the a band and around (tt, 0) in the (3 band. The spin 
structure becomes broader having a spread AQ of the 
size of the pockets, and it still peaks around Q = (0, tt), 
so that Q bridges the two Fermi surfaces. As seen in 
Fig.^(b), the corresponding gap functions have fixed signs 
in each band, and their variation is small. Thus, inter- 
pocket pair scatterings V[k, — k] G A — > V[k', — k'] G B 
(and vice versa) all have significant positive contributions 

tO V ff. 

Correspondingly, as shown in Fig.^(a), the maximum 
eigenvalue A for T = 0.14 starts out around ~ 0.5 for the 
single-band case of t x \ = t X 2 = 1.3, but it increases up to 
~ 0.9 upon decreasing t x \. T c , plotted in Fig.||(b), comes 
close to ~ 0.1 for t x i < 1, which is almost an order of 
magnitude larger compared to typical values for the cases 
with a large connected Fermi surface. 

In order to show that the topology and the size of the 
Fermi surfaces are crucial, we now show results for var- 
ious values of t X 2, fixing t x % — 0.8 (dashed arrow in the 
inset of Figj6|(a)). When t X 2 — t x \ or smaller, the Fermi 
surfaces in both bands are open in the x direction. This 




FIG. 7. Plots of | G a (k) | 2 for 
with t xl = 0.8 and T = 0.14. 



=l.l(a), 1.2(b), and 1.3(c) 



obtained, which, in light of the present results, is be- 
cause the parameter values adopted there do not result 
in small, pocket-like Fermi surfaces. 

In the large U limit, the present Hubbard model should 
tend to the t-J model with alternating t and J. The t- 
J model with alternating J (dimerizpd t-J model) has 
been studied as a spin gapped systemilil There, however, 
the correlation between superconductivity and the topol- 
ogy/size of the Fermi surfaces was not discussed. It is not 
clear at present whether the t-J model with alternation 
only in J has properties similar to those presented here. 
The Hubbard and t-J models on purely ID ladder lat- 
tices have also been studied extensively as spin gapped 
systems. In particular, large enhancement of the pair- 
ing correlation has been obtained when— the-, Fermi level 
lies near the bottom of the upper bandMiill. These re- 
sults seem to be related to the present study, although 
we cannot make a definite conclusion because the FLEX 
approximation cannot be applied to purely ID systems, 
which arc not Fermi liquid but Luttinger liquid. 

At present, it is not clear whether the presence of a spin 
gap is a sufficient or a necessary condition for a super- 
conducting state without nodes of the gap on the Fermi 
surface. Intuitively, however, a gapfull superconductiv- 
ity and a finite spin gap do seem to be consistent with 
each other. This point remains as an appealing future 
problem. 
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FIG. 8. 
T = 0.14. 



A plotted as a function of t X 2 with t x i = 0.8 and 



can be understood by considering the single band situa- 
tion t X 2 — t x i < tyi where the unfolded Fermi surface for 
n ~ 1 is open in the x direction. On the other hand, we 
have seen that there are two small pocket-like Fermi sur- 
faces for t X 2 — 1.3, so each Fermi surface has to change 
its topology between t x i ~ 0.8 and 1.3. In fact, it can be 
seen from FigJ^ that the topology of the Fermi surface in 
the a band changes from a large Fermi surface to a small 
pocket-like Fermi surface at around t X 2 ~ 1.2. A rapidly 
increases corresponding to this change of the topology of 
the Fermi surface, as can be seen from Fig.||. On the 
other hand, A decreases for too large t x2 because the f3 
band lies above the Fermi level in this case.ll3 

Finally, let us discuss the relation between the present 
study and previous Wfs. The Hubbard model on 
coupled-ladder latticesLrH has been studied using FLEX 
approximation. In those studies, T c ~ O(0. Oli) has been 
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